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Abstract 


According to a previous conjecture, spatial and temporal Lyapunov expo¬ 
nents of chaotic extended systems can be obtained from derivatives of a 
suitable function; the entropy potential. The validity and the consequences 
of this hypothesis are explored in detail. The numerical investigation of a 
continuous-time model provides a further confirmation to the existence of the 
entropy potential. Furthermore, it is shown that the knowledge of the en¬ 
tropy potential allows determining also Lyapunov spectra in general reference 
frames where the time-like and space-like axes point along generic directions 
in the space-time plane. Finally, the existence of an entropy potential implies 
that the integrated density of positive exponents (Kolmogorov-Sinai entropy) 
is independent of the chosen reference frame. 
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Lead paragraph 


Out-of-equilibrium systems give rise to a rich variety of phenomena which are often dealt 
with ad hoc methods. Perhaps the main difficulty lying behind the development of general 
tools is the lack of an Ansatz equivalent to the Gibbs-Boltzmann weight, which allows 
estimating a priori the probability of each conhguration of an equilibrium system. Although 
the problem cannot be (easily) overcome, it is tempting to explore the possibility to extend 
and apply concepts like Lyapunov exponents to extract information about the invariant 
measure of a high-dimensional system. For instance, on the one hand it is well known 
that, through the Kaplan-Yorke formula, it is possible to estimate the fractal dimension of 
a hnite (low) dimensional attractor; on the other hand, however, it is unclear how the local 
dynamics arising in different regions of an extended system combine together to determine 
the global invariant measure. It is our opinion that a sensible answer will be given only 
after having extended the Lyapunov analysis, to account for spatial propagation as well 
as temporal instabilities. The chronotopic approach, introduced in Ref. represents a 
tentative construction that goes in this direction, introducing a sort of dispersion relations 
in chaotic systems. In this paper, we test the general validity of the method both with 
numerical simulations and by investigating its internal coherence. 
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I. INTRODUCTION 


Lyapunov exponents, providing information on the evolution in tangent space are very 
useful in that they allow characterizing the invariant measure This is particularly clear 
in the case of hyperbolic systems as revealed by the construction of the Bowen-Ruelle-Sinai 
measure [Q , but is also effectively true for the general and more realistic class of smooth non- 
hyperbolic systems. However, an effective exploitation of these ideas in spatially extended 
systems is hindered by the inhnite-dimensionality of the phase-space and by the existence 
of propagation phenomena that requires a more detailed understanding besides that one 
provided by the knowledge of the usual Lyapunov spectrum [Q. In order to reach a more 
complete comprehension of phenomena occurring in tangent space, two families of Lyapunov 
spectra characterizing temporal, resp. spatial, dynamics of inhnitesimal perturbations have 
been introduced and discussed |]I|. Subsequently P, it has been shown that the two families 
are not independent of one another. More precisely, it has been conjectured that all the 
information on Lyapunov exponents can be obtained from a single observable: the entropy 
potential <h, a function of the temporal and the spatial growth rates A, p, respectively. The 
name “entropy potential” follows from the observation that <I>(A = 0,/i = 0) coincides with 
the density of Kolmogorov-Sinai entropy [§]. However, such a conjecture is based only on 
a few simple examples that can be analytically worked out and on numerical simulations 
performed for a chain of coupled maps. One of the aims of the present paper is to strengthen 
the validity of the conjecture by performing new tests in a more realistic system. 

Moreover, we intend to show a further connection between spatial and temporal Lya¬ 
punov exponents by studying the evolution of perturbations along generic “world-lines” in 
the space-time, i.e. along directions other than the natural space and time axes. The exten¬ 
sion of the usual dehnition of Lyapunov exponents to this more general class of frames, par¬ 
tially discussed in [|^ , is rather appropriate for characterizing patterns with some anisotropy. 
Here, we prove that this seemingly more general class of spatiotemporal exponents can be 
derived from the knowledge of spatial and temporal Lyapunov spectra, which thus conhrm 
to contain all the relevant information. 

The present paper is organized as follows. In Section H we recall the dehnition of the 
entropy potential and present some new test for its existence. Sec. HI is devoted to Lyapunov 
analysis in tilted reference frames i.e. to the dehnition of spatiotemporal exponents, while 
Sec. IV deals with their relationships with more standard dynamical indicators. Some 
conclusive remarks are hnally reported in Sec. V. 


II. LYAPUNOV SPECTRA FROM AN ENTROPY POTENTIAL 

In this Section, we hrst introduce the main classes of models employed in the investigation 
of spatio-temporal chaos. Then, we recall the notion of temporal (TLS) and spatial (SLS) 
Lyapunov spectra with a particular emphasis to their connection with the entropy potential. 
The very existence of the latter is then investigated in the last part, for a specihc continuous¬ 
time model. 

React ion-dihusion systems are among the most relevant models for the study of spa¬ 
tiotemporal chaos. The standard form for the evolution equations is 
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dtu = F(u) + , 


( 1 ) 


where the held u(x, t) is dehned on the domain [0,L], with periodic boundary conditions 
u(0,t) = u{L,t). The nonlinear function F accounts for the local reaction dynamics, while 
the matrix D represents the spatial coupling induced by diffusion. 

A simplihed class of models can be obtained by spatial discretization, i.e. by considering 
a ID lattice of coupled oscillators ^ 

Ui = F(ui) + D (u,+i - 2ui + Ui_i) , (2) 


where the index i labels each site of a lattice of length L (assuming again periodic boundary 
conditions). 

A further simplihcation is achieved by discretizing also the time, i.e by considering cou¬ 
pled map lattices (CML) [^], that in their usual form read as 

<+l = / + (1 “ ’ (3) 


where n is the time index, and e gauges the diffusion strength. The function /, mapping 
a given interval I of the real axis onto itself, simulates the nonlinear reaction process. In 
particular, we will focus on the homogeneous f(x) = rx (mod 1) and logistic f(x) = 
4:x{l-x) CML with < e [0,1]. 

CML models have been also proposed to mimic ID open-flow systems [IT], characterized 
by flux terms in the held equation. For example the model 


Un+l = 


f (^(1 - a )< X (1 - £)< + 


( 4 ) 


with the parameter a bounded between 0 and 1, accounts for the possibility of an asymmetric 
coupling, corresponding to hrst order derivatives in the continuum limit. 

Let us denote with 6{x,t) a generic perturbation that can be assumed to possess an 
exponential prohle both in space and time. 


6{x, t) = a{x, t) exp(—/rx -|- Xt) 


( 5 ) 


Depending whether /r, or A is considered as a free parameter, one can dehne either the 
temporal or the spatial Lyapunov spectrum. Let us hrst discuss the temporal Lyapunov 
spectrum nx{X,fi): it is obtained in the usual way by following the evolution in tangent 
space of a perturbation of the form b{x, t)e^^, where /i is hxed a priori. Obviously, for /r = 0, 
one recovers the standard Lyapunov spectrum. Alternatively, one can consider the temporal 
growth rate A as a free parameter and thereby determine the spatial spectrum n^(A, /i). The 
latter procedure is symmetric to the previous one: fi and A are mutually exchanged, as well 
as X and t, i.e. the tangent dynamics is followed along the spatial direction. 

As a result of the above two approaches one is confronted with a set of four variables: /i, 
A, 77.^, and rix- The main conclusion of is that any two of them can be taken as independent 
variables, the other two variables providing a complete characterization of spatio-temporal 
stability properties. Furthermore, in 0, we have conjectured that both integrated densities 
are not actually independent of one another, but can be obtained from a single function 
$(A,/7), 
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( 6 ) 


= nx 

, 


that is called entropy potential as it coincides with the usual Kolmogorov-Sinai entropy 
density along a suitable line (see Sec. IV). 

Since equivalent descriptions are obtained by choosing any pair of independent variables, 
the appropriate entropy potential is obtained by means of a Legendre transform in the new 
variables. For instance, if one wishes to consider A and nx as new conjugate variables, the 
new entropy potential \h(nA, A) is given by 

^ = Aua - <F , (7) 

and the following relations hold 


dnx 


dfi 


= -n,. 


( 8 ) 


The existence of a potential has been conjectured on the basis of analyticity properties of 
the relationship linking /i and A. The key steps to prove, e.g. (||), are (i) the introduction 
of complex spatial and temporal “growth rates” fl and A, respectively; (ii) the identihca- 
tion (apart from a sign and a multiplicative factor) of Im{X) and with and nx, 

respectively (obviously, Re{\) = A and Re{jl) = p). 

As a matter of fact, we have been able to substantiate the second point only in simple 
cases where the linear stability analysis reduces to the an eigenvalue equation for a fixed 
point or a periodic orbit. In the fully aperiodic regime, there is only numerical evidence 
in a chain of coupled logistic maps, where a test on the existence of a potential has been 
performed by computing a certain circulation integral [Q]. 

In order to test the general validity of our main Ansatz, we consider now a continuous¬ 
time model. Since the generation mechanisms of the multipliers in tangent space are not 
important for our conclusions, we have assumed that they are the result of stochastic pro¬ 
cesses, thus escaping the need to integrate also some nonlinear equations in phase-space. We 
shall refer to a model of coupled oscillators (see Eq. ^), whose evolution in tangent space 
can be written as 


5ui = + D{6ui+i -I- Sui^i) 


( 9 ) 


where the random multipliers mi{t) are indirectly dehned through their power spectrum 


S{u:) = So 


exp 


(a; - ooof 


a 


+ exp 


(a; -F ooof 


a 


( 10 ) 


This choice implies that mi{t) is an analytic function of time and possesses some degree of 
periodicity as one expects to be the case of many deterministic chaotic signals. 

From a numerical point of view, we have decided to compute the derivative appearing in 
(^) through a simple hnite difference (Euler) scheme. This approach, although not the most 
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refined one, allows treating in a consistent way the temporal integration and the recursive 
iteration in space. Indeed, by expressing the perturbation as 6ui(t) = Vi(t)exp {—fii + At), 
Eq. can be rewritten as 


Vi{t + At)e^^* 


Vi{t) + At 


mi{t)vi +D{vi+ie>"+ Vi_ie 


( 11 ) 


which can be easily iterated either in space or time to determine the TSL and the SLS. In 
principle, one should choose a sufficiently small integration time step; in practice the spatial 
recursion becomes soon unfeasible, since the dimension of the “temporal” phase space is 
proportional to 1/At. Accordingly, we have considered a few different values of At, namely 
0.2, 0.1, and 0.05, paying more attention to testing the existence of the potential in each 
case, rather than to computing the actual limit spectra (i.e. to the limit At —> 0). In fact, 
although for a hnite value of At, one cannot assert, rigorously speaking, to have simulated 
Eq. model (|^) can nevertheless be considered as a dynamical system in itself worth 
being studied (actually, it is a CML with a different coupling scheme from that usually 
assumed). 

In Fig. 1^, we have reported the borders of the spectra in the (/i, A)-plane as estimated 
from the SLS and the TLS. The agreement between the two sets of lines for every value 
of At provides a hrst indication that each approximation of the continuous model is in 
itself consistent with the Endings of [jl[. Each border is made of 4 distinct lines: an upper 
and a lower curve representing the set of maxima and minima of TLS; two symmetric lines 
representing the extrema of the SLS. For At —>■ 0, we expect that the first two curves converge 
to an asymptotic shape, while the latter two diverge to ±cxo. In fact, the maximum and the 
minimum spatial divergence rate are limited by the time step At. This is a general feature 
of continuous time systems as already remarked in |]I[. 

The variation of the lateral borders with At indeed confirms our expectations, while the 
dependence of the upper and lower curves (especially at high |/i|-values) indicates, instead, 
that the convergence to the asymptotic form is rather slow: this is clearly confirmed by the 
comparison with the border of the TLS determined for a much smaller value of At (0.00625) 
(see solid curve in Fig. p. 

In any case, the most relevant test that we wanted to perform concerns the the existence 
of an entropy potential <h. This can be numerically done by computing the integral 


nB 

/ {n^,nx)-wds 
A. 


( 12 ) 


along different paths joining generic pairs A, B of points in the (/r. A) plane. If the integral is 
independent of the path, we have a clear indication of the existence of a potential. All tests 
that we have performed for different choices of points and of the time-integration steps have 
been successful!, yielding results that are equal within the numerical accuracy. In Table I, 
we report a subset of the results, all referring to the same end points but different values 
of At. We consider this as a further clear evidence that the potential entropy indeed exists 
and is a rather general feature of at least ID spatio-temporal systems. 


III. SPATIOTEMPORAL EXPONENTS 

In the perspective of a complete characterization of space-time chaos, one should consider 
the possibility of viewing a generic pattern as being generated along directions other than 
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time and space axes. In fact, once a pattern is given, any direction can, a priori, be considered 
as an appropriate “time” axis. Accordingly, qnestions can be addressed abont the statistical 
properties of the pattern when viewed in that way. The extension of the Lyapnnov formalism 
to generic orientations of the space-time coordinates does not follow simply from an abstract 
need of completeness; it also resnlts from the attempt of generalizing the nonlinear time- 
series analysis to patterns. The existence of low-dimensional chaos led many researchers to 
investigate the possibility whether a given irregnlar temporal signal might be the conseqnence 
of a few nonlinearly conpled degrees of freeedom. The existence of space-time chaos leads to 
the eqnivalent qnestion whether a given pattern is the resnlt of a deterministic ID nonlinear 
process. At variance with temporal signals which can be generated only by moving forward 
or backward in time, in the case of patterns, the identihcation of the most appropriate spatial 
and temporal directions is a new and nnavoidable element of the game. For this reason, in 
the next snbsection we shall introdnce the notion of spatiotemporal exponents. 


A. Definitions 

For the moment, we assnme that the pattern is continnons both along space and time 
directions; we shall discnss later how the dehnitions can be extended to CML models. Let 
ns consider a given spatiotemporal conhguration of the held u(x, f), generated, say, by inte¬ 
grating Eq. When arbitrary directions are considered in the {x, t) plane, the coordinates 
mnst be properly scaled in order to force them to have the same dimension. We choose to 
mnltiply the time variable by c, where c is a snitable constant with the dimension of a ve¬ 
locity. Moreover, let'd denote the rotation angle of the tilted frame {x', ct') with respect to 
the initial one {x,ct), adopting the convention that positive angles correspond to clockwise 
rotations. Sometimes it will be more convenient to identify the new frame by referring to 
the velocity v = ctan'd. The limit cases n = 0 ('d = 0) and v = -|-oo ('d = vr/2) cor¬ 
respond to pnrely temporal and pnrely spatial propagations, respectively. The coordinate 
transformation reads as 


ct' = (3 (^ct + -x^ (13) 

x' = {3 {—vt + x) , 


where {3 = 1/y^( 1 -|- n^/c^). The physical meaning of v is transparent: it can be interpreted 
as the velocity in the old frame of a point stationary in the tilted frame (constant x'). 

The new held u(a:', ct') can be thonght of as being the resnlt of the integration of the 
model derived from the original one after the change of variables (|^. Althongh it is not 
obvions whether the invariant measnre in the initial frame is still attracting in the new frame 
(see Ref. [0] for a discnssion of this point), one can anyhow stndy the stability properties by 
linearizing and dehning the Lyapnnov exponents in the nsnal way. 

In CML models, the discreteness of both the space and the time lattice leads to some 
difhcnlties in the practical constrnction of tilted frames. In fact, only rational valnes of the 
velocity v can be realized in hnite lattices (in this case, it is natnral to assnme that the 
lattice spacing is the “same” along the spatial and the temporal directions and, accordingly, 
to set c = 1). Moreover, writing the explicit expression of the model reqnires introdncing 
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different site types. For this reason, we discuss in the following the simplest nontrivial case 
n = 1/2, the generalization to other rational velocities being conceptually straightforward. 

A generic spatial conhguration in the tilted frame is dehned by sites of the spatiotemporal 
lattice (i,n) connected by alternating horizontal (as in the usual case) and diagonal bonds 
(see Fig. By suitably adjusting the relative fraction of the two types of links, all rotations 
between 0 and vr/d can be reproduced. The explicit expression of the updating rule requires 
a proper numbering of the consecutive sites. Moreover, as seen in Fig. it involves the 
“memory” of two previous states. 

Finally, an exact implementation of the mapping rule requires acausal boundary condi¬ 
tions, since the knowledge of future (in the original frame) states is required 0 (this is a 
general problem occurring also in the continuous case). As we are interested in the thermo¬ 
dynamic limit, we bypass the problem by choosing periodic boundary conditions. Such a 
choice has been shown not to affect the bulk properties of the dynamical evolution [0. 

In the updating procedure, two different cases are recognized: the variable u is either 
determined from the past values in the neighbouring sites, or it requires the newly updated 
M-value on the right neighbour (see Fig. For v = 1/2, this can be done by simply 
distinguishing between even and odd sites, 

V"Vi = / (|r»"' + (1 - (14) 

= / (/’ft + (1 - + /’fS/ . 

where i = 1,..., L/2 (L is assumed to be even for simplicity), while 

U+, = XI , (16) 


are additional variables introduced to account for the dependence at time n — 1. Taking into 
account that can be determined from the X and Y variables at time n, the mapping 

can be hnally expressed in the usual synchronous form (X/, F/) —^ (-A/_,_;^, F))_,_^), but with an 
asymmetric spatial coupling with next and next-to-next nearest neighbours. The Lyapunov 
exponents rjj can now be computed with the usual technique . 

In analogy with the original model, we expect again that, in the limit of inhnitely ex¬ 
tended systems, the set of exponents rjjlv) will converge to an asymptotic form. 




(16) 


where n,, = j/L is the corresponding integrated density. We will refer to this function as 
the spatiotemporal Lyapunov spectrum (STLS). In the limit cases u = 0, -|-cx) ("d = 0,7r/2), 
the STLS reduces to the standard temporal and spatial spectrum, respectively. 

The recursive scheme (^) implies an increase of the phase-space dimension by a factor 
(1-I-1/2) (in general 1 -|- v). Actually, as we will argue, these new degrees of freedom are 
not physically relevant. However, for consistency reasons with the original rescaling of the 
spatial variable, we choose to normalize the spatiotemporal density between 0 and 1 -|- u (the 
time units are, instead, left unchanged by the above construction). 

The generalization to asymmetric maps (^ is straightforward: it removes the degeneracy 
V —>■ —V. Numerical results for logistic maps, indicate that the dependence of the positive 
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exponents on the velocity is qnite weak in the fully symmetric case a = 1/2 (for instance, the 
maximum exponent exhibits a 20% variation in the whole v range), while it is remarkable 
for asymmetric couplings. In every case, the negative part of the spectrum sharply changes 
with the velocity. This is consistent with the results obtained for delayed maps in Ref. [0]. 


B. Representation in the (p,A) plane 


Spatiotemporal exponents can be put in relation with /i and A by rewriting the general 
expression for a perturbation in a frame rotated by an angle'd, 

exp(/ix + \t) = exp(pV + \'t') . (17) 

Such an equation induces a rotation of the same angle in the (c/i. A) variables, 

A' = /d (A + vfi) (18) 

/i'=/3 (-(n/c^)A +/i) 

The above equations allow studying the stability with respect to generic perturbations with 
an exponential prohle along x'. For simplicity, we shall consider only uniform perturbations, 

exp(/ix + At) = exp(r 7 T) , (19) 

where the growth rate r] denotes the spatiotemporal exponent. Notice that we have changed 
notations from A' to rj, to understand that the condition /i' = 0 is fulhlled. From the second 
of Eq. (|T8[) , uniform perturbations in the rotated frame correspond to points along the line 
C 


A = c^fx/v 


( 20 ) 


in the (p. A) plane. 

Whenever the evolution of an exponentially localized perturbation of type (|^) is con¬ 
sidered, it is natural to introduce the quantity V = A/p, which can be interpreted as the 
velocity of the front [0. Eq. (^) connects this velocity with that of the rotated frame, 

V = c^jv . (21) 


Therefore, on the basis of dehnition (|T^, {x',cf) can be interpreted as the reference frame 
in which the front associated with the perturbation propagates with an “inhnite” velocity. 
The explicit expression for r] is 

T] = sign(A)^A2 + (cp)2 = X/p . (22) 


Such a relation can be turned into a self-consistent equation for the maximum Lyapunov 
exponent T^max by imposing the constraint that the pair (/i. A) lies on the upper border 
A = Ainax(h) of the spectrum of the Lyapunov exponents, namely 


hmax 



(23) 
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Some ambiguities arise when velocities v > c are considered, since the line C intersects 
Amax(h) in two points as seen in Fig. This phenomenon was already noticed in Ref. 
while discussing the propagation of exponentially localized disturbances in the original ref¬ 
erence frame. Moreover, it has been shown that only the front corresponding to the smaller 
value of /i is stable, except for some cases where a nonlinear mechanism intervenes domi¬ 
nating the propagation process 

At n = c^/14 the two intersections degenerate into a single tangency point. This condition 
dehnes 14, which can be interpreted as the slowest propagation velocity of initially localized 
disturbances [P!^ . 

The extension of Eq. (^) to the rest of the spectrum requires to connect nx and with 
Uri- In the next section, we will show how to perform such a step with the help of the entropy 
potential. Here, we limit ourselves to discuss the structure of the STLS for different values 
of the tilting angle'd. In Fig. we report the borders of the bands, which can be determined 
from the intersections of C with the border &D of the domain of allowed perturbations (see 
Fig. 1 of Ref. and Fig. |^). For ■d = 0 (temporal case) a single band is present but, as soon 
as > 0, a second negative band arises from the intersections with the branch diverging 
to —cxo at p = —Pc- For "d > vr/d, the negative band disappears and a positive band arises 
from the intersections with the branch diverging to -|-cxd with slope v = c. A single band 
spectrum is again recovered for > ^9* = atan(c/14). 

Notice that in symplectic maps, the STLS is symmetric for any value of '& (see Ref. [Q) 
so that positive and negative bands appear and disappear simultaneously. 

It is worthwhile to illustrate some of the above considerations in the simple example of 
the linear diffusion equation 
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dtu = 'yu + Dd^u , (24) 

that arise when dealing with the linear stability analysis of uniform and stationary solutions 
u{x,t) = u* of the scalar version of Eq. (|l|). The expression for Amax(h) can be obtained by 
assuming u{x,t) = exp {fix + At), so that 

Amax(h) = 7 + Dfi^ . (25) 

Accordingly, Eq. (p3D reads as 

^hmax = 1 + D . (26) 

On the other hand, the model equation in the rotated frame can be obtained from the 
substitutions 


dt^ P {dt> - vd^,) (27) 

dx^ + dx^ . 

By introducing the usual Ansatz for the shape of the perturbation, 

u{x\ t') ~ exp [inx + {rf + , (28) 
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separating the real from the imaginary part, and eliminating we obtain the integrated 
density of spatiotemporal exponents 


and the corresponding STLS //(«:, n). Dimensional analysis shows that k is an inverse length, 
as expected for a density of exponents. Notice that Eq. (^) is recovered, by setting k = 
D = 0 in Eq. (|2^). 

In this and more general continuous models, we should remark that the line C intersects 
'^max(h) twice for any arbitrarily small v. This is because the Laplacian operator sets no 
upper limit to the propagation velocity of disturbances. 


IV. FROM THE ENTROPY POTENTIAL TO DYNAMICAL INVARIANTS 

The present section is devoted to establish the consequences of the existence of the 
entropy potential on the Lyapunov spectra and other dynamical indicators such as the 
Kolmogorov-Sinai entropy and the Kaplan-Yorke dimension of the attractor. In order to 
keep the notations as simple as possible, we assume that time and space coordinates are 
scaled in such a way that c = 1. 

A. Spatiotemporal exponents 

The very existence of the entropy potential <I> implies that the Lyapunov spectrum in a 
frame tilted at an angle (recall that "d is the angle from the A-axis) can be obtained by 
computing the derivative of <I> along the direction u = (sin■d, cos'd) in the (/i. A) plane. In 
fact, this is a straightforward generalization of the previous Endings that and nx are the 
derivatives of <h along the /i and A direction, respectively. Accordingly, the STLS is linked 
to the TLS and SLS by the following general equation 

nr,{v,r]) = u= f3[vn^ + nx\ , (30) 

where V = {d^, dx) is the gradient in the {fi, A) plane, and the r.h.s of the above formula is 
evaluated for 


/i = vjSr] (31) 

X = l3v . (32) 

Such a relation can be directly verified for the diffusion equation from Eqs. (8) of Ref. @] and 
Eq. (^). Further, more significative tests have been performed by checking numerically the 
validity of Eq. (|30D in some lattice models involving, e.g., logistic and homogeneous chains, 
the spectra of which are reported in see Fig. ^ (notice that, since the time units have not 
been renormalized in the tilted frame, the factor (3 is no longer needed). 
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B. Entropy 


The connection between the chronotopic approach and the Kolmogorov-Sinai entropy 
Hks was already discussed in [0. Here we recall the main concepts both for the sake 
of completeness and since it can be more effectively phrased with the help of the spatio- 
temporal representation. 

As usual, we shall refer to Pesin formula as a way to estimate the rate of information- 
production from the positive Lyapunov exponents. Although it provides an upper bound, 
there is numerical evidence that the bound is actually saturated in generic models. The 
reference to Pesin’s formula suggests that H^s is an extensive quantity . For this reason, 
it is convenient to introduce the entropy density h\ which, in the thermodynamic limit, can 
be obtained from the integral of the positive part of the Lyapunov spectrum. 

It is natural to extend the definition of entropy to a generic reference frame, 

T^max 

hr,= nr,{'d,7])dr] , (33) 

^ 0 

where the integral is performed along the line C and r/max is the maximum value of t] which 
is reached in the intersection point between C and T). In the limit ■d = 0, the above equation 
reduces to the usual definition of h\, while for "d = 7r/2 it reduces to the “spatial” entropy 

The existence of an entropy potential implies that 




h\ < 'd* 


(34) 


where'd* is the value for which the line C is tangent to V. This is immediately seen by 
combining the observation that follows from an integral along the straight line defining 
the corresponding A-axis with the observation that one of the two extrema is always the 
same (the origin) while the others lie along an equipotential line. 

The independency of of v has an immediate physical interpretation, which we also 
consider as a strong argument in favour of the existence of an entropy potential. The 
Kolmogorov-Sinai entropy density is, in fact, the amount of information necessary to char¬ 
acterize a space-time pattern of temporal duration T and spatial extention S (apart from 
the information flow through the boundaries @]) divided by its area S x T. Therefore, hxs 
is expected to be independent of the way the axes are oriented in the plane, i.e. of the 
velocity n. As a consequence, hri = hxs for all n < 1. 

The above conclusion still holds when the STLS exhibits a positive band as well (which is 
always the case in continuous models), provided that the content of such a band is discarded. 
Accordingly, we can conclude that the new degrees of freedom, associated to the positive 
band, which appear in the rotated frame are just physically irrelevant directions which turn 
the original attractor into a repellor. If n > c^/I4, the two bands merge together and it is not 
anymore possible to distinguish between unstable but irrelevant directions and the unstable 
manifold of the original attractor. Presumably, this means that the repellor is turned into 
a strange repellor with a singular measure along some (all) unstable directions. 
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C. Dimension 


A second important indicator of the “complexity” of a spatiotemporal dynamics is the 
fractal dimension of the underlying measure. An upper bound Dky to it is given by the 
Kaplan-Yorke formula The existence of a limit Lyapunov spectrum, implies that Dky 
is proportional to the system size, so that it is convenient to introduce the dimension den¬ 
sity dxY i, i.e. the number of independent degrees of freedom actively involved in the 
asymptotic evolution per unit length. In the framework of the present paper, it is natural 
to extend the concept of dimension density to a frame oriented in a generic way in the 
space-time plane. A straightforward generalization of the Kaplan-Yorke formula leads to 
the integral equation. 



(35) 


As for the entropy density, Eq. (|35|) can be more easily interpreted with reference to the 
plane. In this representation the entropy potential is 


$ = Xnx + iin^ — d) 


(36) 


In fact, the curve implicitely dehned by the constraint (|35D is the equipotential line C 


*h(n^,nA) = 0 


(37) 


The dimension density dxYiv) can, in turn, be determined from Eq. (^OD at the intersection 
point between C and the image of C in the plane {n^,nx). 

At variance with the entropy density, dxYiv) changes with v (see Fig. P) even if we 
avoid considering the second positive band. In fact, while hrj is an information divided by 
a space-time area, dxYiv) is a number of degrees of freedom divided by a length, measured 
orthogonally to the propagation axis. Thus, at least from a dimensional point of view, it is 
meaningless to compare dxYiv) for different velocities. However, one can reduce temporal 
to spatial lengths by introducing the scaling factor c and, in turn, ask himself how the 
dimension changes with c. It is easily seen that the scaling dependence on c is expressed by 
the following relation. 


dKY{v,Ci)^l+ = dKY{v,C2)^l +, (38) 

The (completely arbitrary) choice of c reflects in different dependences of dxY on n. A 
natural procedure to hx c is by minimizing the dependence of dxY on the observation angle. 
This amounts to choosing the time units in such a way as to make the 2D pattern as isotropic 
as possible. In homogeneous CMLs, the procedure is so effective that a suitable choice of 
c allows removing almost completely the velocity dependence as seen in Fig. where the 
results for the natural value c = 1 are compared with those for c = 3. 

More in general, however, it is not possible to achieve such a complete success. This is, 
for instance, the case of the logistic CML, where the dimension drop for c = 1 is too large 
to be compensated by any choice of c (see Fig. Pd, where the curve for c = 1 is compared 
with the best results obtained for c = -|-cx)). 
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A further indicator which is sometimes useful in characterizing the chaoticity of a given 
extended system is the dimension density of the unstable manifold. This dimension is 
nothing but in the point where ?; = 0, i.e. in the origin, and its expression simply reads 
as 

d„ = /3nA(0,0) . (39) 

being n^(0, 0) = 0. The choice c = +cx) of the scaling factor removes exactly the dependence 
on the orientation of the reference frame. This choice is equivalent to measuring lengths in 
the untilted frame. 


V. CONCLUSIONS AND OPEN PROBLEMS 


In this paper we have provided further numerical evidence for the existence of an entropy 
potential in ID spatially extended dynamical systems. The result appears to follow from the 
possibility to order the Lyapunov vectors according to their average wavenumber that thus 
becomes equivalent to the integrated density of exponents. This interpretation, if conhrmed, 
would imply the possibility to dehne and attribute some meaning to the imaginary part of the 
expansion rates. This line of thought is very reminiscent of the rotation numbers introduced 
by Ruelle in the context of Hamiltonian systems [Q . It is certainly worth to fully explore 
this route in the hope to arrive at a more rigorous justihcation of our theoretical construction 
: work in this direction is in progress. 

Furthermore, we have seen that, among the consequences of the existence of an entropy 
potential, there is the possibility to insert very coherently in this scheme a more general class 
of Lyapunov exponents corresponding to various directions of propagation in the space-time. 
As a result, we have found that the asymmetry between the temporal and the spatial axes can 
be revealed by the dependence of the density of dimension on the orientation of the reference 
frame. This is not true for the Kolmogorov-Sinai entropy density which, as already remarked 
in PI, should be regarded as a “super-invariant” dynamical indicator. 
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FIGURES 
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FIG. 1. Boundaries dV for the model^(P) as obtained from TLS (lines) and SLS (symbols) for 
three different values of the time step At : 0.2 (dotted line and asterisks), 0.1 (dot-dashed line 
and diamonds) and 0.05 (dashed line and circles). The asymprotic boundaries are also reported 
for At = 0.00625 (solid line). 



FIG. 2. 
V = 1/2. 


Lattice implementation of the definition of spatiotemporal Lyapunov exponents for 
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FIG. 3. Plot of the boundary dV and of the line A = U/U in the three velocity regimes for the 
logistic CML £ = 1/3. The three lines refer to the different cases v < 1 (solid), 1 < u < 1/14 
(dot-dashed) and v = 1/14 (dashed). 



0 71/4 nil 


FIG. 4. Boundaries of the STLS versus the tilting angle 'd for the homogeneous chain (' 
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FIG. 5. Comparison between the STLS obtained by direct numerical computation and formula 
(1^ for (a) homogeneous (r = 2) with u = 4/5 and (b) logistic CMLs with u = 3/5 (in both cases 
e = 1/3). 
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TABLES 


Path 

At = 0.2 

At = 0.1 

At = 0.05 

(0,0) ^ (0,1) ^ (2,1) 

-1.1419 

-1.1868 

-1.2303 

(0,0) ^ (2,0) ^ (2,1) 

-1.1358 

-1.1838 

-1.2236 

Relative difference 

0.5 % 

0.2 % 

0.5 % 


TABLE I. Entropy potential for model @ {So = 1/3, a = 0.1, ujq = I, D = 1) computed 
by integrating along two different piecewise linear paths of the (/U, A) plane, and for three different 
time steps. In the last row the relative difference between the two values is reported: relative 
statistical error in the computation of Lyapunov exponents are of the order of 10“^. 
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